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Abstract. This paper describes a general algorithm for finding the commen- 
surator of a non-arithmetic hyperbolic manifold with cusps, and for deciding 
when two such manifolds are commensurable. The method is based on some 
elementary observations regarding horosphere packings and canonical cell de- 
compositions. For example, we use this to find the commensurators of all 
non-arithmetic hyperbolic once-punctured torus bundles over the circle. 

For hyperbolic 3-manifolds, the algorithm has been implemented using 
Goodman's computer program Snap. We use this to determine the commensu- 
rability classes of all cusped hyperbolic 3-manifolds triangulated using at most 
7 ideal tetrahedra, and for the complements of hyperbolic knots and links with 
up to 12 crossings. 



1. Introduction 

Two manifolds or orbifolds M and M' are commensurable if they admit a com- 
mon finite sheeted covering. For hyperboHc n-orbifolds, we can suppose that 
M = H"/r and M' = H"/r', with V and V discrete subgroups of Isom(H"). 
In this paper, we assume that M and M' are of finite volume and of dimension at 
least 3. Then, by Mostow-Prasad Rigidity, commensurability means that we can 
conjugate T by an isometry g such that gTg~^ and T' intersect in a subgroup of 
finite index in both groups. 

Given that the classification of finite volume hyperbolic manifolds up to homeo- 
morphism appears to be hard, it seems sensible to attempt to subdivide the problem 
and start with a classification up to commensurability. Looked at in this way, we 
see a remarkable dichotomy between the arithmetic and non- arithmetic cases. (See 
P5] for the definition of arithmetic hyperbolic manifolds.) 

Define the commensurator of T to be the group 

Comm(r) = {g G Isom(H") | [L : T n gVg''^] < oo}. 

Then T and T' are commensurable if and only if Comm(r) and Comm(r') are con- 
jugate. Geometrically, an element of the normalizer of F in Isom(H") represents a 
symmetry (i.e. isometry) of M = E1I"/F. Similarly, an element of the commensura- 
tor represents an isometry between finite sheeted covers of M; this gives a hidden 
symmetry of M if it is not the lift of an isometry of M (see [27]). 

It follows from deep work of Margulis p4j (see also [35]), that in dimension 
> 3, the commensurator Comm(F) is discrete if and only if F is not arithmetic. 
This means that the commensurability class of a non-arithmetic, cofinite volume, 
discrete group F is particularly simple, consisting only of conjugates of the finite 
index subgroups of Comm(F). In terms of orbifolds, it means that M and M' are 
commensurable if and only if they cover a common quotient orbifold. 

On the other hand, commensurability classes of arithmetic groups are "big" : we 
may well have commensurable F and F' such that the group generated by gTg^^ 
and F' is not discrete for any g. 
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A hyperbolic n-orbifold is cusped if it is non-compact of finite volume. This 
paper describes a practical algorithm for determining when two cusped hyperbolic 
n-manifolds cover a common quotient, and for finding a smallest quotient. For 
non-arithmetic finite volume cusped hyperbolic n-manifolds of dimension n > 3, 
this solves the commensurability problem. 

Section 2 begins with some elementary observations about horoball packings 
and canonical cell decompositions of a cusped hyperbolic manifold. This leads 
to a characterization of the commensurator of a non-arithmetic cusped hyperbolic 
manifold M as the maximal symmetry group of the tilings of H" obtained by lifting 
canonical cell decompositions of M. In Section 3, we use this to determine the 
commensurators of non-arithmetic hyperbolic once-punctured torus bundles over 
the circle. 

Section 4 gives an algorithm for finding the isometry group of a tiling of H" 
arising from a cell decomposition of a hyperbolic manifold, and Sections 5 and 6 
describe methods for finding all possible canonical cell decompositions for a cusped 
hyperbolic manifold. Section 7 contains some observations on commensurability of 
cusps in hyperbolic 3-manifolds which can simplify the search for all canonical cell 
decompositions . 

In 3-dimensions, each orientable hyperbolic orbifold has the form AI — IElI"^/r, 
where F is a discrete subgroup of PSL(2, C) = Isom^(EII^). The invariant trace field 
k(T) C C is the field generated by the traces of the elements of F*^^' = {7^ | 7 S F} 
lifted to SL(2,C). This is a number field if M has finite volume (see [57], [H], 
|22jV The invariant quaternion algebra is the k{T) subalgebra of M2(C) generated 
by F^^^ These are useful and computable commensurability invariants (see |10j . 

ESI). 

For the arithmetic subgroups of Isom(H'^), the invariant quaternion algebra is a 
complete commensurability invariant. In fact for cusped arithmetic hyperbolic 3- 
orbifolds, the invariant trace field is an imaginary quadratic field and the quaternion 
algebra is just the algebra of all 2 x 2 matrices with entries in the invariant trace field 
(see Theorem 3.3.8]); so the invariant trace field is a complete commensurability 
invariant. However most cusped hyperbolic 3-manifolds are non- arithmetic (cf. [B]) 
so other methods are needed to determine commensurability. 

Damian Heard and Oliver Goodman have implemented the algorithms described 
in this paper for non-arithmetic hyperbolic 3-manifolds; these are incorporated in 
the computer program Snap |16| . Using this we have determined the commensura- 
bility classes for all manifolds occurring in the Callahan-Hildebrand- Weeks census 
(|19j. [1]) of cusped hyperbolic manifolds with up to 7 tetrahedra, and for com- 
plements of hyperbolic knots and links up to 12 crossings, supplied by Morwen 
Thistlethwaite (see f20]). These results are discussed in Section 8, while Section 9 
outlines the Dowker-Thistlethwaite notation used to describe links. 

This work has uncovered interesting new examples of commensurable knot and 
link complements (see Examples 2.1 and 2.2 1, and a new example of a knot with 
shape field properly contained in the invariant trace field (see Example 7.11. The 
results have also been used by Button [8] to study fibred and virtually fibred cusped 
hyperbolic 3-manifolds. 

For 1-cusped manifolds we note that "cusp density" (see Section [2]) is a very 
good invariant. We have found only a few examples of incommensurable 1-cusped 
manifolds which are not distinguished by cusp density (see Example 2.3). 

There is also a "dumb" algorithm, based on volume bounds for hyperbolic orb- 
ifolds, which works for any (possibly closed) non- arithmetic hyperbolic 3-orbifold, 
but appears to be quite impractical. If M and M' cover Q with Vol{Q) > C 
then the degrees d, d' of the coverings are bounded hy D — [Vol(A/) /CJ and 
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D' = [Vol(Af')/CJ , respectively. Then if M and M' are commensurable, they 
admit a common covering N of degree at most D' over M and at most D over 
M' . The best current estimate for C for orientable non-arithmetic 3-orbifolds is 
0.041 . . . from recent work of Marshall-Martin [25 . Since Vol(Af) w 2 is typical, 
we would have to find all coverings of M' of degree d! < 50. This means finding 
all conjugacy classes of transitive representations of 7ri(il/') into 6*50, a group with 
around lO^"* elements! 

Acknowledgements: We thank Ian Agol for pointing out a simplification to our 
method of determining the commensurability of Euclidean tori, Gaven Martin for 
information on current volume bounds, and Walter Neumann for several interesting 
discussions on this work. We also thank Alan Reid, Genevieve Walsh, and the 
referee for their helpful comments on the paper. 

2. The Commensurability Criterion 

We use the following terminology throughout this paper. A set of disjoint 
horoballs in H" is called a horohall packing, and a cusp neighbourhood in a hy- 
perbolic 71-orbifold is one which lifts to such a horoball packing. 

Lemma 2.1. The symmetry group of a horohall packing m H" is discrete whenever 
the totally geodesic suhspace spanned by their ideal points has dimension at least 
n — 1. 

Proof. Let {gi} be a sequence of symmetries of the packing converging to the iden- 
tity. Choose horoballs i?i , . . . , i?„ whose ideal points span a totally geodesic sub- 
space H of dimension n—\. For i sufficiently large, we can assume that gi{Bi.) — 
for k = 1, . . . ,n. But this implies that these gi fix H pointwise. Since the only such 
isometries are the identity, and reflection in H, the sequence must be eventually 
constant. □ 

Lemma 2.2. Let M = IHI"/r be a finite volume cusped hyperbolic orbifold. The set 
of parabolic fixed points ofT spans H". 

Proof. The set of parabolic fixed points is dense in the limit set of T which equals 
the whole of the sphere at infinity. □ 

Lemma 2.3. Let M , M' be finite volume cusped hyperbolic orbifolds. Then M and 
M' cover a common orbifold Q if and only if they admit choices of cusp neighbour- 
hoods lifting to isometric horoball packings. 

Proof. If M and M' cover Q, choose cusp neighbourhoods in Q and lift to M and 
M' . These all lift to the same horoball packing in H", namely the horoball packing 
determined by our choice of cusp neighbourhoods in Q. Conversely, both M and 
M' cover the quotient of H" by the group of symmetries of their common horoball 
packing which, by Lemmas |2 . 1 1 and |2 .2| is discrete. □ 

We can define the cusp density of a 1-cusped hyperbolic orbifold M as follows. 
Since M has only one cusp it has a unique maximal (embedded) cusp neighbourhood 
U. The cusp density of M is Vol{U) /Vol{M). Since the cusp density of any orbifold 
covered by M is the same, it is a commensurability invariant of orbifolds with 
discrete commensurator. 

Choosing a full set of disjoint cusp neighbourhoods in a non-compact finite vol- 
ume hyperbolic n-manifold M determines a "Ford spine." This is the cell complex 
given by the set of points in M equidistant from the cusp neighbourhoods in two 
or more directions. Cells of dimension n — k contain points equidistant from the 
cusp neighbourhoods in fc + 1 independent directions {k = 1, . . . , n). This spine can 
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also be seen intuitively as the "bumping locus" of the cusp neighbourhoods: blow 
up the cusp neighbourhoods until they press against each other and flatten. 

Dual to the Ford spine is a decomposition of M into ideal polytopes, generically 
simplices. The ideal cell dual to a given 0-cell of the Ford spine lifts to the con- 
vex hull in H" of the set of ideal points determined by the equidistant directions. 
We call the cell decompositions that arise in this way canonicaZ|^ For a 1-cusped 
manifold the canonical cell decomposition is unique. It is shown in ^ that a finite 
volume hyperbolic manifold with multiple cusps admits finitely many canonical cell 
decompositions . 

Theorem 2.4. Cusped hyperbolic n-manifolds M and M' cover a common orbifold 
if and only if they admit canonical ideal cell decompositions lifting to isometric 
tilings o/H". 

Proof. If M and AI' cover Q, choose cusp neighbourhoods in Q and lift them to 
M, M' and H". Constructing the Ford spine and cell decomposition in H" clearly 
yields the lifts of those entities from both M and Af corresponding to our choice 
of cusp neighbourhoods. 

Conversely, observe that the symmetry group of the common tiling gives an 
orbifold which is a quotient of both manifolds. □ 

Remark 2.5. We can omit the word 'canonical' in the above theorem. The proof 
is unchanged. 

The previous theorem gives the following characterization of the commensurator. 

Theorem 2.6. Let M — IHI"/r be a finite volume cusped hyperbolic n-manifold 
with discrete commensurator. Then Comm(r) is the maximal symmetry group of 
the tilings o/H" obtained by lifting canonical cell decompositions of M; it contains 
all such symmetry groups. 

For manifolds with discrete commensurator, we can now define a truly canon- 
ical ideal cell decomposition as follows. Find Comm(r) as in the above theorem. 
Choose equal volume cusp neighbourhoods in ElI"/Comm(r). Lift them to M and 
take the resulting canonical cell decomposition of M. Two such manifolds are com- 
mensurable if and only if their truly canonical cell decompositions give isometric 
tilings of H". 

Choosing maximal cusp neighbourhoods in IEII"/Comm(r) also gives a canonical 
version of cusp density for multi-cusped manifolds. 

Theorem |2.6| is the basis for the algorithms described in this paper. Canoni- 
cal cell decompositions can be computed by the algorithms of Weeks described in 
[33j and implemented in SnapPea |34j . In Section 4 below we give an algorithm 
for finding the isometry groups of the corresponding tilings of H". Combining 
this with Theorem |2.6| gives an algorithm for finding commensurators of 1-cusped 
non-arithmetic hyperbolic n-manifolds. In Sections 5 and 6 we extend this algo- 
rithm to multi-cusped manifolds, by describing methods for finding all canonical 
cell decompositions. 

For hyperbolic 3-manifolds, these algorithms have been implemented by Heard 
and Goodman. (These are incorporated in "find commensurator" and related 
commands in the program Snap [16j). We conclude this section with some examples 
discovered during this work. 



The term is not really ideal since, for manifolds with multiple cusps, there are generally 
multiple canonical cell decompositions depending on the choice of cusp neighbourhoods. 
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2.1. Example: a 5-link chain and friends. The following five links have com- 
mensurable complements, as shown using our computer program. In the first three 
cases at least it is possible to 'see' this commensurability. 




The first of these is the 5-link chain C5. Thurston [SH Chapter 6] explains how to 
obtain a fundamental region for the hyperbolic fc-link chain complements: we span 
each link of the chain by a disk in the obvious manner. The complement of the 
union of these five disks is then a solid torus. Once the link is deleted, the disks and 
their arcs of intersection divide the boundary of the solid torus into ideal squares 
A, B, C,D,E as shown, with cusps labelled a, 6, c, c?, e. 



(6) (e 



a) (b) (c 



m 



CD 



O D 



n 



— @ — ® — © — ® — ^ 

A T bTcTdTeT/ 

— — — — ^ 



A hyperbolic structure is given by taking two regular pentagonal drums with 
ideal vertices and adjusting their heights to obtain (ideal) square faces. Glue two 
drums together as shown, identify the top with the bottom via a ^ rotation, and 
glue faces as indicated. Edges are then identified in 4's, two horizontal with two 
vertical. It is easy to check that the sum of dihedral angles around each edge is 27r 
so this gives a hyperbolic structure, since the angle sum is tt at each ideal vertex of 
a drum. 
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It is clear from the symmetry of the picture that the drums are cells in the canon- 
ical cell decomposition obtained by choosing equal area cusp cross sections. We 
remark that Neumann-Reid show that this link complement is non-arithmetic in 
EZl Section 5]. 

Now change the 5-link chain by cutting along the shown disk, applying a half 
twist, and re-gluing to obtain our second link. In the complement, this surgery 
introduces a half turn into the gluing between the A-faces. Edges are still identified 
in 4's, two horizontal with two vertical, but there are now only four cusps. 




If we repeat the process on the second disk shown we obtain the third link. Again 
this corresponds to changing the gluing pattern on our two drums. 

Since these link complements are non-arithmetic, the tiling of H'^ by pentagonal 
drums covers some canonical cell decomposition of each one. Since their volumes 
are the same, each one decomposes into two pentagonal drums. It should therefore 
be possible, in each case, to find 5 ideal squares meeting at order 4 edges, cutting 
the complement into one or two solid tori. We leave this as a challenge for the 
reader. 

2.2. Example: commensurable knot complements. Commensurable knot com- 
plements seem to be rather rare. Previously known examples include the Rubinstein- 
Aitchison dodecahedral knots |3] and examples due to a construction of Gonzales- 
Acufia and Whitten [T^ giving knot complements covering other knot complements. 
For example, the —2, 3, 7 pretzel knot has 18/1 and 19/1 surgeries giving lens spaces. 
Taking the universal covers of these lens spaces gives new hyperbolic knots in 
whose complements are 18- and 19-fold cylic covers of the (— 2, 3, 7)-pretzel com- 
plement . 

Our program finds a pair of knots "9n6" and "12n642" , having 9 and 12 crossings 
respectively, whose complements are commensurable with volumes in the ratio 3 : 4. 
Walter Neumann has pointed out that these knots belong to a very pretty family 
of knots: take a band of k repeats of a trefoil with the ends given m half twists 
before putting them together. E.g. {k,m) — (3,2): 
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The half-twists are put in so as to undo some of the crossings of the trefoils (allowing 
the above projection to be rearranged so as to have 9 crossings). The pair of knots 
found by our program correspond to (fc, m) = (3, 2) and (4, 1). 

To see that these knots have commensurable complements we find a common 
quotient orbifold. In each case this is the quotient of the knot complement by its 
symmetry group; these are dihedral groups of order 12 and order 16 respectively. 

The picture of 9n6 above shows an obvious axis of 2-fold symmetry; below left is 
the quotient, which is the complement of a knot in the orbifold with singular set 
an unknot labelled 2. By pulling the knot straight, we see that this is an orbifold 
whose underlying space is a solid torus with knotted singular locus. 



If we arrange the singular locus on a torus parallel to the boundary of the solid 
torus we see 3 clasps, 3 strands in the (vertical) core direction, and a strand with 
slope 2/1. The view from inside the solid torus looking towards the boundary is 
shown below. (For the knot (4, 1) we would see 4 clasps, 4 strands in the core 
direction and a strand with slope 1/1.) 



The solid torus with its singular set has three 2- fold symmetries whose axes intersect 
the solid torus in 6 arcs, each passing perpendicularly through the core like a skewer, 
with symmetry group dihedral of order 6. The ends of the arcs are shown as stars 
above. 

The quotient orbifold is obtained by taking a shce of the solid torus between 
two axes and folding closed the top and bottom disks like books. The result is a 
ball with the axes giving two unknotted arcs of order 2 in the singular set, running 
out to the boundary (which is now a (2, 2, 2, 2)-pillowcase orbifold). The original 
singular set gives an arc linking the other two, so that the whole singular locus is 
an 'H' graph labelled with 2's. 

The last three pictures show what happens to the singular locus in one slice of 
the solid torus as we fold. We begin with the annulus in the bottom l/6th of the 
previous figure, redrawn after twisting the bottom. This bounds a solid cylinder 
with the singular locus as shown in the middle figure. Folding along the top and 
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bottom (and expanding the region slightly) gives the final result. 




We leave it to the reader to draw similar pictures for the knot with {k,m) — (4, 1) 
and verify that the result is indeed the same orbifold. Alternatively, Orb [TS] or 
SnapPea [34, can be used to verify that the appropriate dihedral covers of the final 
orbifold give the complements of the knots with (fc, m) — (3, 2) and (4, 1). 

We remark that Walter Neumann has found an infinite family of new examples 
of pairs of commensurable knot complements in the 3-sphere; this example is the 
simplest case. 

2.3. Example: cusp horoball pictures. Figure [l] shows the horoball packings of 
two 1-cusped census manifolds ml37 and ml38 as seen from the cusp. Using Snap 
|16j . we find that the commensurability classes of these two equal- volume mani- 
folds are indistinguishable by cusp density or invariant trace field. Their maximal 
horoball packings and canonical cell decompositions are however different. (For 
example, the edges joining degree 4 vertices in the following cusp diagrams are all 
parallel for ml37, but not for ml38.) 




Figure 1. The maximal horoball packing and canonical cell de- 
composition as seen from the cusps of manifolds inl37 and ml38. 
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3. Example: punctured torus bundles 

The bundles over with a once-punctured torus as fibre provide an interesting 
family of 1-cuspcd hyperbolic 3-manifolds with known canonical triangulations. 
By analysing how symmetries of the lifted triangulations of the universal cover 
appear when viewed from the cusp we obtain (in all the non- arithmetic cases) 
strong constraints on what symmetries may be possible. We then show that they 
all come from symmetries of the manifold. This leads to the following: 

Theorem 3.1. Let M — M^/G be an orientable non- arithmetic hyperbolic 3- 
manifold which is a once-punctured torus bundle over S^. Then M has no "hidden 
symmetries", i.e. the commensurator of G is the normalizer of G in Isom(IHI^). 

Let F denote a once-punctured torus, and let : f ^ f be an orientation 
preserving homeomorphism. Let 

1 ^ X [0, 11 

(x,0) - (v9(a;),l) 

be the mapping torus of tp. Identifying F with (M^ — Z^)/Z^, we have that (p is 
isotopic to an element of SL(2,Z); since M depends only on the isotopy class of 
tp we assume tp G SL(2,Z). Then M is hyperbolic whenever p is hyperbolic, i.e. 
when tp has distinct real eigenvalues; and M^pi are homeomorphic if and only 
if p and p' are conjugate. 
Define matrices 

( 1 1 )' ^= ( 1 

For each word w in the symbols L,R define p^ G SL(2,Z) as the corresponding 
matrix product. 

Lemma 3.2. Each p G 5'L(2,Z) is conjugate to ±pm for some word w in the 
symbols L, R. The sign is unique and w is determined up to cyclic permutations of 
its letters. 

Let M — where p ~ dzp^u is hyperbolic. The so-called monodromy trian- 
gulation T of M has one tetrahedron for each letter in w and gluings determined 
by w and the sign. It is nicely described in [12] and [14]. It follows from work of 
Lackenby [5T] that T is the canonical ideal cell decomposition of M. Other proofs of 
this result have recently been given by Gueritaud ([2], [T^) and Akiyoshi, Sakuma, 
Wada and Yamashita (see 5 ). 

The intersection of T with a (small) torus cross section of the cusp of M, lifted to 
its universal cover M^, gives the (lifted) cusp triangulation Tq of M. Note that edges 
and vertices of Tq correspond to edges of T seen transversely or end-on respectively. 
SnapPea [53] provides pictures of these cusp triangulations: see Figure [3] for an 
example. 

We need two things: the first is a combinatorial description of Tq in terms of w; 
the second is an understanding of which edges and vertices of To correspond to the 
same edges of T in M. Both are outlined briefly here: for detailed explanations we 
refer the reader to the Appendix of [12] and Sections 3 and 4 of [I^ . 



3.1. The monodromy triangulation T. The triangulation T is built up in layers 
by gluing tetrahedra according to the letters of w. We begin with an almost flat 
ideal tetrahedron projecting onto a punctured torus; the tetrahedron has edges 
a, b, c, c_ identified as shown below. 



10 



OLIVER GOODMAN, DAMIAN HEARD, AND CRAIG HODGSON 




For each successive letter L or i? in the word w, we attach a tetrahedron to the 
top of the previous tetrahedron, as shown below in the cover (M^ — Z^) x M. 





After using all the letters of w, the final triangulation of the fibre F differs from 
the initial triangulation by the monodromy tp, and we can glue the top and bottom 
together to obtain an ideal triangulation T of M. 



3.2. Combinatorial description of Tq. Now consider the induced triangulation 
of a cusp linking torus (i.e. cusp cross section) in M. Each tetrahedron contributes 
a chain of 4 triangles going once around the cusp as shown in Figure [2j 




Figure 2. The chain of triangles around a cusp coming from one tetrahedron. 



In the triangulation To of this lifts to an infinite chain of triangles forming a 
(vertical) saw-tooth pattern. 
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Each chain is glued to the next in one of two ways depending on whether the 
letter is an L or an R. 




Before gluing we adjust slightly the "front" triangles of the first chain and the 
"back" triangles of the second, so as to create two horizontal edges. After stacking 
these chains of triangles together we have a decomposition of Tq into horizontal 
strips which can be described combinatorially as follows. 

Make a horizontal strip out of triangles corresponding to the letters of w. for 
each L add a triangle with a side on the bottom edge of the strip and a vertex 
at the top; vice-versa for i?; repeat infinitely in both directions. For example, for 
LRRLR we have: 




Fill the plane by reflecting repeatedly in the top and bottom edges of the strip. 
Then the cusp linking torus is the quotient of the plane by the group Go generated 
by vertical translation by four strips, and horizontal translation by one period of 
the strip (composed with an extra vertical translation by two strips if 1^9 = —'^^w)- 




Figure 3. The cusp triangulation of M = M^p^ where w = 
LRRLR. The triangulation of AI has five edges, labelled 0-4, which 
appear as both edges and vertices of the cusp triangulation. 
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3.3. Edge/vertex correspondence in Tq. To understand which edges of Tq are 
identified in T, we define a "direction" on each horizontal strip: right to left on the 
first strip, left to right on the strips below and above it, and so on, so that adjacent 
strips have opposite directions. Then we have: 

Lemma 3.3. Each edge e of Tq which crosses a strip bounds two triangles, one of 
which, S say, lies on the side of e given by the direction on the strip. Then e and 
the opposite vertex of S give the same edge ofT. Horizontal edges of Tq correspond 
with the edges at the opposite vertices of both adjacent triangles. 

This result is illustrated in Figure |4] and will be important in the arguments 
below. 



a 


b 


c 






c\ 


'i 




4 








\' 
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b 


c 


a' 





reflected strip 
initial strip 



Figure 4. Edge identifications in strips of Tq. 



Proof. The boundary between two tetrahedra in T is a punctured torus consisting 
of two ideal triangles, homotopic to a fibre of A/. In Tq the boundaries give edge 
cycles, lifts of a cycle of 6 edges going once around the cusp in the vertical direction 
as in Figure [2j Here are two cycles: 



Cycles meeting a given strip form cither "L's and F's" or "backwards L's and F's." 
This gives the directions: L's and F's read left to right! 

We see from Figure [2] that the edges of a cycle are homotopic to three edges 
of T, a, h, c say, cyclically repeated. The vertices are the same three edges of T 
arranged so that if an edge is a, its vertices are h and c and so on (as in Figure 
|4]). Thus for any (forwards or backwards) L or F: the edge at the down-stroke (a 
strip-crossing edge) equals the edge seen at the other end of its horizontal stroke. 
Further, the edge seen at any horizontal stroke equals the edge at the other end of 
the down-stroke. This proves the lemma. □ 

3.4. Reduction to Tq. Let G be the image of the holonomy representation of 
7ri(Af) in Isom(H'^) so that M = /G as usual. Let a: be a fixed point of some 
maximal parabolic subgroup of G which we identify with the group of translations 
Go defined above. We regard symmetries and hidden symmetries as isometrics of 
H^. Since M is 1-cusped we can compose any symmetry or hidden symmetry with 
an clement of G to obtain an isometry which fixes x. Thus if Af has a symmetry. 
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it is represented by a symmetry of Tq; if it has a non-trivial hidden symmetry, it is 
represented by a symmetry of Tq which does not come from a symmetry of M . 

Since our combinatorial picture of Tq is not metrically exact , we only know that 
symmetries and hidden symmetries act as simplicial homeomorphisms (S.H.) of Tg. 

Lemma 3.4. Simplicial homeomorphisms oJTq preserve horizontal strips whenever 
w 7^ (LR)"^ or {LLRR)"^ as a cyclic word, for any m > 0. 

Since vertex orders in Tq are all even we can define a straight line to be a path 
in the 1-skeleton which enters and leaves each vertex along opposite edges. A strip 
consists of a part of Tq between two (infinite, disjoint) straight lines such that every 
interior edge crosses from one side of the strip to the other. 

Proof. If an edge of Tq hes on the edge of a strip, not necessarily horizontal, the 
opposite vertex of the triangle which crosses that strip must have order at least 6: 
the straight line going through this vertex has at least the two edges of the triangle 
on one side. 

If there exists an S.H. which is not horizontal strip preserving, every vertex of Tq 
will lie on a non-horizontal edge which is an edge of some strip (namely the image 
of a previously horizontal strip edge) . 

Suppose there is a vertex of order 10 or more, corresponding to 3 or more adjacent 
Us or _R's in w: 



Here none of the non-horizontal edges shown can lie on a strip edge because they are 
all opposite vertices of order 4. So in this case all S.H.'s must preserve horizontal 
strips. 

Suppose there is a vertex of order 8, corresponding to LL or RR in w: 



Now the only non-horizontal edges shown which can be strip edges are the vertical 
ones. Vertex orders on this vertical edge alternate 4,8,4,8. So if an S.H. maps a 
horizontal strip edge to this vertical one, vertex orders along some horizontal edge 
must be 4,8,4,8, . . .. This determines Tq and hence w as {LLRR)"^ . (The reader 
can verify that this particular Tq admits a tt/2 rotation.) 

Finally suppose there is no vertex of order > 6. It follows immediately that w is 
{LR)™ and Tq is the tiling of the plane by equilateral triangles. (Again, this admits 
S.H.'s which are not horizontal strip preserving.) □ 

Lemma 3.5. Simplicial homeomorphisms ofTa coming from symmetries ofT pre- 
serve the horizontal strip directions whenever w ^ (LR)"^ or (LLRR)"^ . 

Proof. An S.H. of Tq which comes from a symmetry of T obviously preserves the 
order of every edge of T. Thus if we label each edge of Tq with the order of the 
corresponding edge of T, the labels must be preserved. 

The order of an edge of T corresponding to a vertex of Tq is simply its order 



as a vertex of Tq. By Lemma 3.3 we can use the strip directions to label the 
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corresponding edges of Tq. Given a vertex of order 10 or more we have: 




The direction on a horizontal strip cannot be reversed because this would swap an 
edge of order 4 with one of order 6 or more. 

Suppose now the maximum vertex order is 8. For each order 8 vertex we have: 



If there is an S.H. which reverses strip direction then the other two sides of this 
diamond figure, wherever it appears, must also be labelled with 4's. 



The 4 on the upper right edge implies that the right hand vertex of the diamond 
has order > 8, hence 8. This gives another diamond figure adjacent to the first one. 
The argument can be repeated giving the pattern associated with w — {LLRRY^ . 
If there are no vertices of order > 6, w = (LR)"^. □ 



Proof of Theorem 3.1 It is shown in [7 that the only arithmetic orientable hy- 
perbolic punctured torus bundles are those for which w — LR, LLR (or LRR) 
or LLRR, or powers of these, having invariant trace fields Q(V — 3), Q(-\ /— 7) and 



3.4 



and 



3.5 



show 



Q('\/— 1) respectively. Thus, in the non-arithmetic cases. Lemmas 
that all symmetries of Tq which come from hidden symmetries of M are represented 
by simplicial homeomorphism which preserve the strips and strip-directions. 

Note that there are two "sister" manifolds for each w depending on whether 
the monodromy is (p^ or —ipw] the triangulations of however are the same, 
depending only on w. 

It is now easy to see that the possible symmetries of Tq , modulo the translations 
in Go, are restricted to the types listed below. We show that each, if it occurs at 
all, comes from an actual symmetry of M . 

(1) Shifting up or down by two strips 
^ ' 1 from F x±^p. 



-1 
is central in SL(2,Z).) 



This is realized by the symmetry 
(Note that ' ^ 



S'l to itself 







-1 
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(2) If w is a power, w = m™, then Tq admits a horizontal translation by £{u) 
triangles, where i{u) is the length of the word u. In this case 

and 1 X Tn^im gives a symmetry of order m (where re denotes rotation of 
the unit circle by 6*). For the —^p^ case we replace one of the (/j„'s by 
—^u- After rotating we also have to apply — 1 x 1 to -F x [0, 1/m]. 

(3) If w is palindromic, i.e. w and its reverse w' are the same as cyclic words, 
then To admits a rotation by tt about a point on one of the strip edges. This 

1 




-1 



1 X refljS^) 



X refl{S^), where refl{S^) denotes a reflection of 




X 1 




Note that conjugation by 



takes L to L ^ and R to R ^. 



1 

-1 

(4) If £{w) is even and rotating w by a half turn swaps L's and R's, Tq has a 
glide reflection mapping a strip to itself, exchanging the top and bottom 

of the strip. This is realized by ^ ^ ^ ^ x 1 since conjugation by this 

matrix swaps L and R. More explicitly, let w = uv where v is u with i's 
and i?'s interchanged. The symmetry is 





followed by r,r on the factor. This symmetry is orientation reversing. 

(5) If £{w) is even and reversing w swaps Us and i?'s (wc might say that w 
is anti-palindromic) then Tq admits a glide reflection with vertical axis, 

shifting everything up by one strip. This is realized by 1 



refl{S^) since conjugation by 



takes L to R ^ and R to L ^. 



1 
-1 

This symmetry is orientation reversing. We leave the picture as an exercise 
for the reader. 

□ 

Clearly if M admits any two of the symmetries 3-5 it admits the third which is 
a product of the other two. Thus if w is not a power, M may have one of 5 possible 
symmetry groups. 
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4. Algorithm for finding isometries of tilings 

To apply our commensurability criterion, Theorem |2.4| we need a way to deter- 
mine when tihngs of H", arising from ideal cell decompositions, are isometric. We 
pursue this question in a slightly more general setting. 

Let M and M' be geometric n-manifolds with universal cover X" — E" or 
H", each with a given finite decomposition into convex polyhedra. We wish to 
determine whether the tilings T and T' which cover these two cell decompositions 
are isometric. (The elements of T (resp. T') are convex polyhedra in X" which 
project to polyhedra in the decomposition of M (resp. M').) 

Necessary and sufficient conditions for the tilings to be isometric are as follows. 

(1) There is an isometry i of which maps an element of T isometrically 
onto an element of T' . 

(2) Whenever i maps P ^ T isometrically onto P' € T' , and F is a (codimen- 
sion 1) face of P, i maps the neighbour of P at isometrically onto the 
neighbour of P' at i{F). 

Sufficiency follows from the fact that we can proceed from any tile in T to any other 
by a finite sequence of steps between neighbouring tiles. 

Let S and S' denote the polyhedra decomposing M and M' respectively. 

Let 6 denote the set of all triples {j,p,p'), where p & S, p' £ S' , and j is an 
isometry carrying p onto p'. We can find in a finite number of steps. Condition 1 
above is equivalent to Q being non-empty. 

We say that {j,p,p') G is induced by an isometry i of X" if there exist P E T 
projecting to p and P' £ T' projecting to p' such that i carries P isometrically 
onto P' and the restriction induces j. Let / be a face of p and let q and q' be the 
neighbours of p and p' at / and j (/) respectively. The restriction of j to / induces 
an isometry between certain faces of q and q'. If this extends to {jq,q,q') we say 
that {j^PtP') extends across f to {jq,q,q')- If {j,P,p') is induced by i then {jq,q,q') 
will be also. 

Condition 2 is equivalent to the following: whenever {j,p,p') is induced by i, 
and / is a face of p, then {j,p,p') extends across /. 

Theorem 4.1. With the above notation, the tilings T and T' are isometric if and 
only ij there exists a non-empty subset I of & such that every element of I extends 
across each of its faces to yield another element of I. 

Proof. If T and T' are isometric with isometry i, simply let / be the set of elements 
of O induced by i. 

Conversely, suppose we have ICQ having the stated properties. Choose any 
element {j,p,p') £ I and any isometry i of X" which induces it. Clearly condition 1 
above is satisfied. For condition 2 note that if i maps P onto P' and this induces 
{j^P:P') G I then, because this extends across all of its faces, i maps neighbours of 
P isometrically onto neighbours of P' and these too induce elements of /. Therefore 
condition 2 is satisfied with all the induced triples belonging to /. □ 

An algorithm for finding such a subset / C or establishing that none exists 
is straightforward. If Q is empty, stop; there is no such subset. Otherwise choose 
any element of Q and put it in /. For each element of / check if we can extend 
across all faces. If we can't, remove all elements of / from Q and start again. If 
we can, and every element we reach is already in /, stop and output /. If we can 
extend across faces of all elements of /, and we obtain new elements not yet in /, 
add those new elements to I and check again. 

Each set / found by the algorithm represents an equivalence class of isometries 
carrying T onto T': choose an isometry i inducing any element of /; then isometries 
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are equivalent if i' = g'ig where g is a, covering transformation of M and g' is 
a covering transformation of M'. 

To find the symmetry group of a tiling T we apply the algorithm with T = T' . 
M = M' etc. Let F denote the group of covering transformations of M. Repeating 
the algorithm until 6 is exhausted we obtain finitely many equivalence classes of 
symmetries. These represent the double cosets of T in Symm(T). Together with F 
they generate Symm(r). 

Remark 4.2. Each set I found by this algorithm also represents a common covering 
N of M and M' constructed as follows: For each element r = {j,p,p') G / let Pr 
denote a copy of p. Let N be the disjoint union of the polyhedrapr as r varies over 
I, with the following face identifications. Each face f of p appears as a face fr of 
Pr- Whenever r extends across f to s = {jq,q,q'), identify face fr of Pr with face 
fs ofqs. 

4.1. Example. Let M be the Euclidean torus obtained by gluing a 1 by 3 rectangle 
along opposite pairs of edges. Let M' be another such torus obtained by gluing 
a 2 by 1 rectangle. Subdivide each into unit squares as shown, so that S, the 
subdivision of M, equals {A, B, C} while S' = {P, Q}. Both S and S' lift to tiUngs 
of the plane by unit squares. 



p 




6 consists of 48 elements since there are 8 ways of isometrically mapping one 

Tinit square onto another and 6 combinations of squares to be mapped. Suppose 
we start with [tAP, A, P) where Iap denotes the translation carrying A onto P 
in the picture. We can extend {tAP,A,P) across the marked edge (whose image 
in Af is also shown) to obtain {tsp, B , P). Continuing until we have extended 
across every available edge we obtain 6 elements of Q which we can abbreviate to 
{tAP,tBP,tcp,tAQ,tBQ,tcQ}- These give rise to the following common covering 
of M and M'. 



t 

AP 


t 

BP 


t 

CP 


t 

AQ 


1 

BQ 


t 

CQ 



Suppose instead we begin with a clockwise rotation tap through 90° carrying A 
onto P. Then we obtain a different common cover. 



r 




r 




r 




AP 


BQ 


CP 


AQ 


BP 


CQ 
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Suppose finally we subdivide one of the squares in S" into two triangles. Now the 
tilings are clearly not isometric. 





































/ 














P 








A 


B 


C 










/ 





















ABC 




The algorithm might start with Iap but when it tries to extend this across the 
marked edge it will find that the induced mapping from an edge of A to an edge of 
Q does not extend to any isometry from A to Q. 

4.2. Combinatorial construction of coverings. Mostow-Prasad rigidity often 
allows us to work in a purely combinatorial setting; this motivates the follow- 
ing. The algorithm described above can be viewed as a method for constructing 
a common cover of two manifolds which admit suitable subdivisions into polyhc- 
dra. The algorithm works equally well for topological spaces M and M' admitting 
n-dimensional polyhedral decompositions (i.e. obtained by gluing n-dimensional, 
Euclidean or finite volume hyperbolic, convex polyhedra pairwise along their (n — 1)- 
dimensional faces, such that open faces of dimension q < n embed). We then let 
each element {j,p,p') G 6 represent a combinatorial isomorphism between p and 
p' . (We can barycentrically subdivide the decompositions of M and M' in order 
to realize the j's as piecewise linear homeomorphisms in a canonical way.) Now, if 
the algorithm is successful it constructs, in general, a branched covering of M and 
M' , branched over the codimension 2 skeleta of M and M' . 

When M and M' are PL-manifolds, we would like to construct an unbranched 
cover which is also a manifold. It is sufficient to check that the cover we construct 
does not branch over any codimension 2 cell of M or M' . For if this is the case 
we proceed by induction on q to show that our cover is not branched over any 
codimension q cell of M (or M'). By construction, there is no branching over 
cells of codimension and 1 , and we assume that there is no branching over cells of 
codimension 2. So suppose g > 3 and there is no branching over cells of codimension 
< q—1. Then the link of a codimension q cell in the cover is a connected unbranched 
covering of the link of a codimension q cell of M. But since the latter is a (q — 1)- 
sphere, hence simply connected since g > 3, the covering is a homeomorphism. 

In the case of most interest to us, namely when M and M' are cusped hyperbolic 
3-manifolds, Mostow-Prasad rigidity ensures that if TV topologically covers both M 
and M' then the induced hyperbolic structures on N are isometric. 

Frequently, an ideal cell decomposition of a hyperbolic 3-manifold consists en- 
tirely of ideal tetrahedra. If this is the case for M and M' , then our algorithm 
will always find a common cover, branching over the edge sets of the two mani- 
folds (elements of 6 will always extend over their faces because all tetrahedra are 
combinatorially equivalent) . This will indicate that the two manifolds are commen- 
surable only if the covering is in fact unbranched. When an element of O induces a 
mapping of an edge e of (the ideal cell decomposition of) M onto an edge e' of M', 
the resulting cover will be unbranched if and only if the order of e (the number of 
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tetrahedra to which glue around it) is equal to the order of e'. For example, this 
has the following corollary. 

Theorem 4.3. // a finite volume hyperbolic S-manifold M admits a decomposition 
into ideal tetrahedra such that the order of every edge is 6, then M is commensurable 
with the complement of the figure-^ knot. 

5. Enumerating canonical cell decompositions I 

For 1-cusped hyperbolic manifolds with discrete commensurator we now have 
all the ingredients of an effective method for testing commensurability. For multi- 
cusped manifolds we need a way to search through the (finite) set of all canonical 
ideal cell decompositions. In this section and the next we describe two alternative 
approaches to this problem, while in Section [7] we show how the search can be 
restricted for greater efficiency. 

The first approach is very simple minded: if we can bound the degree with which 
M covers its commensurator quotient, we can enumerate all cusp cross sections 
which could possibly cover equal area cross sections in the quotient. Such degree 
bounds can be obtained from estimates on the the minimum volume of cusped 
non- arithmetic hyperbolic 3-manifolds ([2S], [I], [2H])- 

The second approach is geometric and may be of more theoretical interest since 
it truly finds aZZthe canonical cell decompositions. In fact it associates with a given 
c-cusped M, a convex polytope in whose fc-dimensional faces, for < A: < c, are 
in 1-1 correspondence with canonical cell decompositions of M. 

Let M be a hyperbolic orbifold with (horoball) cusp neighbourhoods Ci , . . . , Cm 
and let iV be a degree d quotient of M with corresponding cusp neighbourhoods 
Ci, . . . , c„. Let TT : M ^ TV be the covering projection and let Cj(i) — ■n{Ci). If we 
choose horospherical cross sections in N with "area" (i.e. codimension 1 volume) 
equal to 1, the area of dCi is equal to the degree with which Ci covers Cj(i). The 
sum of the areas of the Ci covering Cj will be d for each Cj . 

Thus, in order to find a possible degree d quotient of M via the canonical cell 
decomposition, we can enumerate the possible integer area vectors as follows. For 
each n, 1 < n < m, and each partition of {1, . . . ,m} into n non-empty subsets 
/i, . . . , J„, enumerate all area vectors (ai, . . . , a^) such that each is a positive 
integer and X^iG/ ~ -^'-"^ j = ^, ■ ■ ■ 

Example: if ttI = 3 we have partitions {{1,2,3}}, {{1,2}, {3}}, {{1,3}, {2}}, 
{{1},{2,3}} and {{1}, {2}, {3}}. The first partition admits \{d - \){d - 2) area 
vectors, the next three admit d — 1 each, and the last, just one. 

We can enumerate area vectors corresponding to possible quotients of degree d 
as follows. Any quotient of M of degree d has n < m cusps and has a corresponding 
area vector (ai, . . . ,am) with sum X]"=i = if we fix an integer D > m 

and enumerate all positive integer area vectors summing to at most D, we will 
find all quotients of degree d with n < m cusps such that nd < D. In particular 
this will include all quotients with at most to — 1 cusps provided d < D/(m — 1). 
Since the canonical cell decomposition is determined by the ratio of the areas, any 
TO cusped quotient is found using the area vector (1, . . . , 1) with ai — m < D. 
So we will find all canonical cell decompositions arising from quotients of degree 
d<D/{m^l). 

6. Enumerating canonical cell decompositions II 

Let M be a hyperbolic n-manifold with c > cusps. Then the set of possible 
choices of (not necessarily disjoint) horospherical cross sections dual to the cusps of 
M is parametrized by the vector of their areas (i.e. (n — l)-dimensional volumes) 
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in — {{vi, . . . , Wc) G K'^ I wi > 0, . . . , Wc > 0}. In fact it turns out to be more 
convenient to parametrize by the (n — l)th root of area, a quantity we will refer to 
as size. Multiplying a size vector v e M^g t>y a constant A > has the effect of 
shifting the corresponding horospherical cross sections a distance log(A) down (i.e. 
away from) the cusps. 

For each v € M^p "^"^ obtain a Ford spine. It is clear from the definition that if 
we choose a set of disjoint horospherical cross sections and then shift them all up 
or down by the same amount we get the same Ford spine. Thus the set of possible 
Ford spines is parametrized by the set of rays in MS,q, or equivalently, by points in 
the open (c — l)-dimensional simplex S — {{vi, . . . , Vc) G IR>o | wi + • • • + Wc = 1}- 
For a 1-cusped manifold the Ford spine is unique. 

Dual to each Ford spine is a canonical cell decomposition D{v). As we vary 
V S the decomposition changes only when the combinatorics of the spine changes. 
To better understand this dependence we now review an alternative approach to 
defining D{v), namely the original one of Epstein-Penner in [IT| . 

We work in Minkowski space E"^^ with the inner product * defined by 

x*y = xiyi + . . . + x„j/„ - x„+i?/„+i 

for X — (xi, . . . , Xn+i), y = • • ■ , Vn+i) in IR"^"'^. Then hyperbolic space H" is 
the upper sheet of the hyperboloid x ^ x = —1, and the horospheres in H" are 
represented by the intersections of hyperplanes, having light-like normal vectors, 
with H". Each such hyperplane H has a unique Minkowski normal n such that 
X G H if and only if a; * n = — 1 . 

Let r denote the group of covering transformations of H" over M. Each size 
vector V gives rise to a F- invariant set of horospheres in H". The resulting set 
of normals in Minkowski space is invariant under the action of the group F. The 
convex hull of this set of points, which we shall refer to as the Epstein-Penner 
convex hull, intersects every ray based at the origin passing through a point in the 
upper sheet of the hyperboloid. The boundary of this convex set is a union of 
closed convex n-dimensional polytopes having coplanar light-like vertices. Epstein 
and Penner llj show that these project to a locally finite, F-invariant set of ideal 
polyhedra in H", which in turn project to a finite set D{v) of ideal hyperbolic 
polyhedra in M . 

Starting with a given ideal hyperbolic cell decomposition D of M and a size 
vector V, we proceed next to describe necessary and sufficient conditions for D to 
be the canonical cell decomposition D{v), as in [33j and [30j. 

Let C be a cell of D. Then v determines a horospherical cross section to each 
ideal vertex of C. We lift C and this choice of horospheres to H". In Minkowski 
space, this gives a convex (Euclidean) n-dimensional polytope whose vertices are 
the hyperplane normals for these horospheres. Whenever C is not a simplex, it is 
necessary to add the condition that these vertices are coplanar. If this is satisfied 
for all the non-simplicial cells of Z?, we can lift each cell to an n-dimensional poly- 
tope in Minkowski space with vertices corresponding to the choice of horospheres 
determined by the size vector v. li C and C" are neighbouring cells in D it is nec- 
essary that the angle between neighbouring lifts into Minkowski space be convex 
upwards. Together these conditions are also sufficient to imply D = D{v). 

These conditions can be expressed as a set of linear equations and linear inequal- 
ities on the entries of the size vector v. 

Proposition 6.1. Let D be an ideal hyperbolic cell decomposition of a cusped hy- 
perbolic n-manifold M with c cusps. Then there exist matrices L]j and Fd with c 
columns such that, for v £ IR>o; D(v) ~ D if and only if L^v ~ and Fov > 0. 
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Note: here v is written as a column vector, and the condition Fov > means that 
each entry of the vector F^v is positive. 

Proof. Let Vi denote the entry of v corresponding to the ith cusp of AI. Let iij be 
the vertex representative for the jth vertex of C hfted to Minkowski space for the 
choice of horospheres given by v = (1, . . . , 1). Then for an arbitrary size vector v 
the corresponding representative is rij/ucQ), where c{j) is the cusp of the jth vertex 
of C. 

The coplanarity condition on a non-simphcial cell C gives a set of linear equa- 
tions satisfied by v, one for each vertex of C in excess of rt + 1, as follows. Let 
Ny be a Euclidean normal to the hyperplane containing {no/wc(o): • ■ • : '^n/i'c(n)} 
such that {nj/vc{j)) ■ — 1 for j — 0,...,n, where • denotes the Euclidean 
dot product. Writing Mc for the inverse of the matrix with rows we obtain 
Ny = M(7(wj,(Q), . . . , which is a linear function of v. For j > n, Y\.j/v^(^j^ 

belongs to this hyperplane if and only if Wj ■ Ny — — 0, which is linear in v. 
The full set of constraints for C gives a matrix equation Lev ~ 0. 

The convexity condition at an {n — l)-cell f of D, being the common face of 
n-cells C and C , can be expressed as follows. Let Ny, as above, be the defining 
normal for the hyperplane containing the lift of C determined by v. Let n'f./vc(k) 
be a vertex of an adjacent lift of C", not in the lift of /. This vertex lies above the 
hyperplane if and only if (n'j,/wc(fc)) • iVi, > 1, or equivalently n'^, • Ny — Vc(k) > 0- 
We refer to the left-hand side of this inequality as the tilt d± f oi v and express the 
condition as FfV > 0, where Ff is a suitable row-vector. (Note that the sign of our 
tilt function is opposite to that of [33] and [3T].) 

Finally, concatenate the matrices Lq into a matrix L^j and the rows Ff into a 
matrix F^. □ 

Let Vd denote the set of u e MS,q such that L^v = and Fjjv > 0. We call 
this the parameter cell of D since it contains all cusp size parameters v such that 
D{v) = D. Each v E MS,g belongs to a parameter cell, namely V{v) = 'Po{v)- The 
parameter cell Vd is non-empty if and only if is a canonical cell decomposition. 

It is shown in [Ij that the number of canonical cell decompositions is finite. 
Therefore MS,g is a union of finitely many parameter cells. 

Proposition 6.2. Each v € MS^q can be perturbed to obtain a nearby vector v' such 
that 'P{v') has dimension c and V{v) is a face ofV{v') (or equals V{v) if this has 
dimension c). 

Proof. Let D = D(y). If Ld is zero (or empty) then V{v) is an open subset, hence 
a c-dimensional cell and we just set v' = v. 

Otherwise, perturb v such that it leaves the linear subspace determined by L^v — 
0. For a small perturbation the Epstein-Penner convex hull changes as follows: no 
dihedral angle between adjacent n- faces goes to tt but some non-simplicial n- faces 
may be subdivided if their vertices become non-coplanar. It follows that may 
lose rows and F^ may gain rows. Let D' be the new decomposition. Then Ld' ^ Lu 
because L^'v' = while L£iv' ^ 0. Repeat until L^i is zero (or empty). Then 
Viy') has dimension c. 

Now V belongs to a face of Vd' since it satisfies FfV > for each face / common 
to D and D', and Ffiv = for each face /' of D' not in D. The former condition 
amounts to F^v > 0. We have to show that the latter is equivalent to L^v — 0. 
But that is equivalent to the coplanarity of the lifted vertices of each n-cell of D. 
Such a cell may be subdivided by new faces /' in D' . Then the vertices will be 
coplanar at v if and only if the tilt at each subdividing face is zero, i.e. if and only 
if FfV = for the subdividing faces. □ 
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The above proposition implies that each face of a parameter cell is another 
parameter cell; the decomposition corresponding to a parameter cell is a refinement 
of the decompositions corresponding to its faces. 

Remark 6.3. It is tempting to suppose that the canonical cell decomposition of 
M corresponding to a c dimensional parameter cell must consist entirely of ideal 
simplices hut this need not he the case. In general we can have cell decompositions 
with non-simplicial cells such that L]j is a zero matrix. (For example, this occurs 
for the Borromean rings complement — see Section \6. 1\ helow.) 



We now have the following algorithm for finding all canonical cell decompositions. 
First we find all the c-dimcnsional parameter cells. 

(1) Choose an arbitrary v E ]R>o- 



(2) Perturb v if necessary, as in the proof of Proposition 6.2 so that V{v) has 
dimension c, and add it to our list of cells. 

(3) If the closure of the cells we have found so far does not contain the whole of 
MS^Q, choose a new v not in the closure of any cell found so far and repeat 
step 2. 

By the finiteness result quoted above, this algorithm eventually terminates. We 
can then enumerate all canonical cell decompositions by enumerating the faces of 
all dimensions of the cells V{v). 

While the computational geometry involved in implementing the above algorithm 
is certainly possible, it is not particularly nice. We explain a refinement which gives 
a little more insight and an algorithm which is easier to implement. 

For a decomposition D of M, let S^i denote the row vector obtained by adding 
together the rows of Fo . We define the tilt polytope of M to be the set oi v E ^>o 
such that ■ v < 1 for all canonical cell decompositions D of M. 

Proposition 6.4. The tilt polytope T of M is hounded. The parameter cells Vd 
of M are the cones over the origin of those faces of T which are not contained in 

Proof. We show that the closure of a c-dimensional parameter cell V d has bounded 
intersection with T. Let u be a unit vector in Vd- Then since v is not contained 
in every face of Vd: Sd • v > 0. The length of any multiple of v contained in T 
is bounded by ■ v). Since this is continuous in v, and the set of such v is 

compact, Vd is bounded. Since T is a union of finitely many such sets it is 
bounded. 

Let us write Hd for the half-space S E'^ | S^j ■ x < 1}. Then T is the 
intersection of all the Hd's with M^q. We will show that: if v helongs to a c- 
dimensional parameter cell Vd, and Vd' is any other parameter cell, then the ray 
generated hy v leaves Hd hefore it leaves Hd' . 

It will then follow that a ray in Vd penetrates the (non-empty) face of T gener- 
ated by Hd. Since a ray not inVo belongs to the closure of some other parameter 
cell, it does not leave Hd first and therefore does not pass through the same face 
of T. Since cones on the lower dimensional faces of a (c — 1) -dimensional face 
of T are the faces of a c-dimensional parameter cell, the result then follows from 
Proposition |6.2| 

It remains to show that Hd cuts off any ray in Vd closer to the origin than Hd', 
for all parameter cells Vd' ^ Vd- Equivalently, for v CzVd, ^d ■ v > ^d' ■ v. 

Firstly, let Vd and Vd' be any two parameter cells such that Vd' is a face of 
Vd, and let v belong to Vd- The rows of Fd' are a proper subset of the rows of 
Fd, and since Ff - v > for each row, T,d ■ v > 'Ed' ■ v. If instead Vd is a face 
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of Vd', then the rows of Fu' omitted from Fd are precisely those for which Ff ■ v 
vanishes. Therefore in that case S^i ■ v = Y,o' ■ v. 

Next, let Vd be a c-dimensional parameter cell, and let Vd' be arbitrary, with 
D' ^ D. Choose v € Vd and v' £ Vd'- Let Vd^ ■ ■ ■ ^'Pd^ be the parameter cells 
through which the straight line Vt (1 — t)v + tv' passes for < t < 1 (so that 
Di = D and Dm — D'). For vt in VDi, ■ Vt > "^Di+i ■ Vt while for vt € VDi+i, 
Sd- ■ Wt < T.Di^-1 ■ Vt- Since the difference between these terms is (afRne) linear in t, 
the former inequality must hold for all lesser values of t, in particular, when vt = v. 
Note also that the first such inequality is strict, namely, I^Di > - v. It follows 
that TiD ■ V > ^d' ■ V for arbitrary D' . □ 

Let To be a polytope resulting from the intersection of ES,q with some of the 
half-spaces Hd defined in the above proof. If Tq D T, some face A = Tq O OHd 
of To will contain a point v G M>q not in T, and thus not in Tn OHd, nor in the 
cone on this, Vd- Therefore Ldv or FfV < for some row Ff of Fd- This 
gives a test for when Tq properly contains T; when satisfied, it yields a new half- 
space Hd{v) whose intersection with Tq is strictly smaller. After a finite number of 
intersections we arrive at Tq = T . See Figure [5] 




Figure 5. Tq is a partially computed tilt polytope. The face of 
To containing v G Vd has a vertex w not itiVd- Therefore D{w) 
gives another face of T. 



The computational geometry involved in the above is relatively straightforward. 
By using homogeneous coordinates we can treat an unbounded region, such as K>o, 
as a polytope with some vertices "at infinity". The face A of To, as defined above, 
is the convex hull of those vertices v of Tg satisfying S^i ■ v = \. If any of these 
satisfy Ldv ^ or FfV < for some row Ff of Fd we conclude that To ^ T. If 
such V lies in i9K>o perturb it a httle to bring it inside R^q before determining 
a new half-space Hd(v)- 

6.1. Example: The Borromean rings complement. Let M be the comple- 
ment of the Borromean rings in 
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Then M is an arithmetic hyperbohc 3-manifold with 3 cusps. It may be realized 
by gluing two regular ideal hyperbolic octahedra in the following pattern (see [32]). 




Letters indicate face identifications; cusps are numbered 0,1,2. Equal area cusp 
cross sections give rise to symmetrically placed ideal vertex cross sections in the 
two octahedra. It follows that the corresponding ideal cell decomposition consists 
of precisely these two octahedra. We will first compute its parameter cell. 

Let us position the tiling of H'^ by ideal octahedra such that one of them, call it 
C, has a lift with vertices at {(±1,0,0), (0, ±1,0), (0,0, ±1)} in the projective ball 
model. When the cusp cross sections all have equal area and contain the centre of 
the octahedron, their Minkowski normals are np = (1, 0, 0, 1), rii = (0, 1, 0, 1), = 
(0,0,1,1), n3 = (-l,0,0, 1), n4 = (0,-1,0,1), ng = (0, 0, -1, 1). A neighbouring lift 
of the other octahedron, call it C", has vertex representatives {ng, ni, n2, rig, 114, n'g} 
where is the Minkowski metric reflection of in the hyperplane spanned by 
{no,ni,n2}. Then nj, = (1, 2, 2, 3), n^^ = (2, 1, 2, 3), n[; = (2, 2, 1, 3). (See the left 
hand side of Figure [6j) 



[1:0:0] 




Figure 6. Initial cell decomposition of M into ideal octahedra 
and the corresponding projective parameter cell. 



For size vector v — {vq, vi, W2)* the resulting ideal vertex representatives become 
Ui/vi, Ui^^/vi and n^_^_^/vi for i = 0, 1, 2. The hyperplane containing no/wo, ni/vi, 
11.2/ V2 and n^/vQ is {x | x • iV„ = 1} where is calculated as in the proof of 
Proposition |6.1| giving 
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Then 114 /vi and ni^/v2 lie in this hyperplane if and only if ■ Ny 
(2, -2, 0)v and ■ ~ V2 = ^ (2, 0, -2)v. Hence ic = ' ^ 




-2 







2 - '-^y 

symmetry, Lc' is the same. Next we compute the tilt of /o, the face between C 
and C . This will be positive if n'^/i'a lies above the hyperplane defined by N^, i.e. 
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if n'^- Ny-vo = (-2,2,2)w > 0. Thus Ff^ = (-2,2,2), and the tilt is positive at 
w = (1, 1, 1)*. By symmetry, the other faces of the octahedra also have positive tilt. 
The parameter cell of this decomposition {v \ vq ~ vi = V2} is, projectively, a point 
(as shown in the right half of Figure [6]). 

Suppose next we increase vq slightly. The faces of C still have positive tilt. We 
claim that C is subdivided into 4 tetrahedra containing the edge no,n3 as shown 
in the left of Figure [7] below. 




Figure 7. Decomposition into ideal tetrahedra and corresponding 
parameter cell. 



Let /i be the triangle with vertices {no,n2,n3}. Then /i has positive tilt if 
and only if n4 ■ Ny — vi = (2, —2, 0)v = Ff\v > 0. Letting /2 be the triangle with 
vertices {no,ni,n3}, we obtain Ff^ — (2,0,-2). By symmetry, the tilts of the 
other two faces shown above are the same. Four more triangles, having the same 
tilts, subdivide C . Therefore M is subdivided into 8 simplices in the parameter 
cell shown in the right of Figure [Tj 

Increasing vq, eventually Ff^v = 0, so that the faces of C and C" vanish from 
the decomposition. The four simplices around the edge ni,n2 become an octahe- 
dron. Further increasing vq so Ff^v < 0, we claim that this octahedron splits at 
the square dual to the edge ni , 112 into two square-based pyramids with vertices 
{no, iIq, 113, na, rii} and {no, Hq, rig, n3, as shown in the left of Figure [s] 




Figure 8. Decomposition into square-based pyramids and corre- 
sponding parameter cell. 



Let define the hyperplane containing {no/uo, ni/ui, ng/uo, ng/wo}- Then 



( 



\ 






/ 



V. 



Let Ci be the square-based pyramid with vertices {no, ni, ns, ng, nf,}, where no = 
( — 1, 2, 2, 3) is the Minkowski reflection of njj in the plane containing ni, n2, n4, ns. 
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We find that Lc^ = (0,0,0). Let fz be the triangle with vertices {nQ^n^.n'^}. 
Then F/, = (1, —1, —1) and Ff.^v > 0. Therefore M is divided into 4 square-based 
pyramids in the parameter cell shown in the right of Figure [8] as mentioned in 
Remark 16.31 




Figure 9. All parameter cells for canonical decompositions of the 
Borromean rings complement 

By symmetry, the full set of parameter cells is projectively as shown in Figure 
|9] On the three line segments containing the centre point, each of the original 
octahedra is subdivided into two square based pyramids. 

7. COMMENSURABILITY OF CUSPS 

The number of canonical cell decompositions that can arise for a multi-cusped 
manifold can be quite large, placing practical limits on the usefulness of our meth- 
ods. We show next how it is often possible to greatly reduce the complexity of 
computing the commensurator of a manifold. 

A horospherical cross section of a cusp in an orientable hyperbolic 3-manifold 
is a Euclidean torus, well defined up to similarity. We can position and scale a 
fundamental parallelogram in C such that one vertex lies at the origin and the 
two adjacent edges end at 1 and a point z in the upper half-plane. Such a z is 
called a cusp shape parameter. Alternative choices of fundamental parallelogram 
yield parameters differing by the action of SL2Z by Mobius transformations. By 
choosing a suitable fundamental domain for the action of SL2Z on the upper half- 
plane, we can make a canonical choice of shape parameter for each cusp. 

If one cusp covers another, there is an induced covering of Euclidean tori. It 
follows that their cusp shape parameters are related by the action of an element of 
GL2Q. Let us call cusp shapes commensurable if they are so related. As we shall 
see shortly, it is not hard to determine when two cusp shapes are commensurable. 

Suppose we are trying to determine the (discrete) commensurator quotient Q 
of Af. If M has cusps of incommensurable shape, these necessarily cover distinct 
cusps of Q. Therefore any assignment of cusp neighbourhoods in Q will yield a Ford 
spine and tiling whose symmetry group is the whole commensurator (the group of 
covering transformations of Q). It follows that we can start by choosing arbitrary 
horospheres in each of a set of representatives for the commensurability classes of 
cusps of M. Then as we vary our choices of horosphere in the remaining cusps we 
will be sure to find a tiling whose symmetry group is the commensurator of M. 
The easiest case is when no two cusps of M are commensurable. Then any choice 
of horospheres at all will do. 

The harder case is when M has multiple commensurable cusps. If the symmetry 
group of M is non-trivial, it may act by non-trivially permuting some of these 
cusps. Conceptually we should first divide M by its symmetry group and then find 
the commensurator of this quotient. For practical purposes it is hard to work with 
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non-manifold quotients. Instead, whenever two cusps are related by a symmetry 
of M, we choose symmetrically equivalent cusp neighbourhoods. Let c be the 
number of orbits under the action of Symm(Af) on the cusps. Let d be the number 
of distinct commensurability classes of cusp shape. Then the parameter space of 
relative horosphere positions we need to search, in order to find a tiling whose 
symmetry group equals the commensurator, has dimension c — d. 

Our algorithm is really only a slight modification of the algorithm of Weeks [33] 
for finding the symmetries of a cusped hyperbolic manifold M = H'^/F. In order to 
find the symmetries we need only consider the canonical cell decomposition arising 
from a choice of cusp neighbourhoods such that all boundary tori have equal area. 
Then the symmetries of M are the symmetries of the lifted tiling that normalize F. 
Equivalently, they are the symmetries of the tiling for which the covering described 
in Remark |4.2| has degree one. 

Returning to the question of when two cusp shapes are commensurable, we note 
first that cusp shapes of M belong to the invariant trace field of M. But if k is any 
number field, and a, a' are irrational elements of k, they are related by an element 
of GL2Q if and only if 

(1) {ca + d)a' = aa + b 
is soluble for a, 6, c, d e Q such that 

(2) ad~bc^ 0. 

We can replace ^ with the condition that a, b, c, d are not all zero, since for 
a, a' ^ Q, ^ follows automatically from (jlj and the fact that {l,a} are lin- 
early independent over Q. Regarding fc as a finite dimensional vector space over 
Q we see that ([ij has non trivial solutions if and only if {1, a, a' , aa'} are linearly 
dependent over Q. In particular, if [fc : Q] < 4 all irrationals are commensurable in 
this sense. We thank Ian Agol for pointing out this condition. 

7.1. Example: a knot with cusp field not equal to invariant trace field. 

An interesting example uncovered during this work is the complement of the knot 
12n706 shown below. This has one torus cusp with shape parameter z — 6i generat- 
ing a cusp field Q(i) which is strictly contained in its invariant trace field q{i,V3). 
This answers a question of Neumann-Reid in [27 , who asked whether the figure 
eight knot and the two dodecahedral knots of Aitchison-Rubinstein were the 
only such exampletF] 




8. Experimental Results 

We have implemented the algorithms described here and used them to com- 
pute the commensurability classes of all 4929 manifolds in the Hildebrand- Weeks 
census of cusped hyperbolic 3-manifolds [5] and all 7969 complements of the 8614 



Alan Rcid informs us that Nathan Dunfield has found another example; the 15 crossing knot 
15nl32539. 
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hyperbolic knots and links up to 12 crossings. After replacing non-orientable man- 
ifolds (in the 5-census) by their double covers and removing duplicates we obtained 
a total of 12783 orientable manifolds. Tables of all the results are available at 
|http://www.ms .unim elb . edu . au/~snap/ 

In the process it was necessary to identify the arithmetic manifolds and arrange 
them into their own separate commensurability classes. Criteria for arithmeticity 
and for commensurability of arithmetic manifolds are given in [lOj and [23j . In fact, 
a cusped manifold is arithmetic if and only if its invariant trace field is imaginary 
quadratic and it has integer traces. Cusped arithmetic manifolds are commensu- 
rable if and only if they have the same invariant trace field. They are therefore 
classified by the discriminant d of the invariant trace field Q{y/d). 

There were 142 arithmetic manifolds in six commensurability classes. Table [l] 
lists some of the arithmetic manifolds found with discriminant —3, —4 or —7, and 
all with discriminant —8,-11 or —15. 



d manifolds 

"~ mOOO, m002, m003, m004=4al, m025, m203=6a5=lln318, sll8, 

s961, 8a39, 10a280, 10nl30, 10nl43, 10nl55, lln539, 12a3285, ... 
"~ mOOl, ml24=8nl0=10nl39, ml25, ml26, ml27, ml28, ml29, s859, 

vl858, 8n8=9n34=10nll2, 9n36, 10a242, 10n74, lln545, . . . 
~ m009, mOlO, s772, s773, s774, s775, s776=6a8=8n7, 8a25, 8a37, 

10n73=12n968, 10nll3, lln498, 12n3068, 12n3078, 12n3093, ... 
~~ v2787, v2788, v2789, 9a73, 9a74, 12a3292, 12a3296, 12n2625= 

12n2630, 12n2972, 12n3088, 12n3098=12n3099. 
"^11 12a2126, 12a2961, 12a3039, 12a3230, 12a3295. 
-15 12a3169, 12a3273, 12a3284, 12a3300, 12a3307, 12a3308. 

Table 1. Selected arithmetic manifolds. 



The naming of manifolds in the tables is as follows. Manifolds whose names 
begin with m, s or v belong to the 5, 6 or 7 tetrahedra census of cusped manifolds 
respectively. The rest are knot and link complements in the form <numher of 
crossings> <alternating or non-alternating> <index in table>. The link tables 
were provided by Morwen Thistlethwaite and are included with current versions of 
Snap [ig and Tube |I7]. 

The remaining 12641 manifolds were non-arithmetic, falling into 11709 commen- 
surability classes. A few of them are shown in Table [2] 

Column headings are as follows: s.g. is the order of the symmetry group; c.d. 
is the degree of the manifold over its commensurator quotient; c.vol is the volume 
of the commensurator quotient, optionally followed by c.cl, commensurability class 
numbered from 0, when incommensurable manifolds are listed with the same com- 
mensurator volume. Thus manifolds are commensurable if and only if they have 
the same entry in this column. The invariant trace field is described by its degree, 
discriminant and a number specifying which root of the minimum polynomial gen- 
erates it (with sign corresponding to choice of complex conjugate); nc,cnc gives the 
number of cusps in the manifold and the number of cusps in the commensurator 
quotient; cusp density is computed using equal area cusps in the commensurator 
quotient. 

The first group of manifolds have the smallest commensurator volume found 
(among non-arithmetic manifolds) and are the link complements which appeared 
in Section 2.1 They all have 'hidden symmetries,' i.e. commensurabilities not aris- 
ing from the symmetry group of the manifold. The total number of non-arithmetic 
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manifold 


s.g. 


c.d. 


c.voljC.cl 


field 


nc,cnc 


cusp density 


10a297 


20 


40 


0.365076519 


4,1025(2) 


5,1 


0.642619992 


10a291 


4 


40 


0.365076519 


4,1025(2) 


4,1 


0.642619992 


10a277 


4 


40 


0.365076519 


4,1025(2) 


3,1 


0.642619992 


12n2492 


4 


40 


0.365076519 


4,1025(2) 


3,1 


0.642619992 


12n2899 


4 


40 


0.365076519 


4,1025(2) 


4,1 


0.642619992 


12nll89 


2 


2 


7.175483613,0 


7,-76154488(-2) 


2,2 


0.589830477 


12nll90 


2 


2 


7.175483613,1 


7,-76154488(-2) 


2,2 


0.589830477 


12nl481 


2 


2 


7.175483613,2 


7,-76154488(-2) 


2,2 


0.589830477 


12n2348 


2 


2 


7.175483613,3 


7,-76154488(-2) 


3,3 


0.644747497 


12n2580 


2 


2 


7.175483613,4 


7,-76154488(-2) 


3,3 


0.631898787 


m045 


4 


4 


0.818967911 


3,-107(-2) 


1,1 


0.608307263 


m046 


4 


4 


0.818967911 


3,-107(-2) 


1,1 


0.608307263 


v3379 


8 


8 


0.818967911 


3,-107(-2) 


2,1 


0.608307263 


v3383 


4 


8 


0.818967911 


3,-107(-2) 


3,1 


0.608307263 


v3384 


8 


8 


0.818967911 


3,-107(-2) 


2,1 


0.608307263 


12al743 


8 


16 


0.818967911 


3,-107(-2) 


2,1 


0.608307263 


v3376 


4 


4 


1.637935822,0 


3,-107(-2) 


2,1 


0.690189995 


v3377 


4 


4 


1.637935822,0 


3,-107(2) 


1,1 


0.690189995 


v3378 


4 


4 


1.637935822,0 


3,-107(2) 


1,1 


0.690189995 


12a2937 


8 


8 


1.637935822,1 


6,-1225043(1) 


3,2 


0.608307263 


9a94 


4 


24 


0.575553268 


4,144(1) 


3,2 


0.844133714 



Table 2. Selected non-arithmetic manifolds. 



manifolds having hidden symmetries was 148. The next group of manifolds shows 
that incommensurable manifolds are not always distinguished by cusp density. The 
third group of manifolds includes manifolds whose classes are distinguished by in- 
variant trace field but not by cusp density, and manifolds with the same commen- 
surator quotient volume but different invariant trace fields. One should not get the 
impression that cusp density is a poor invariant: in fact, among the 11278 cusp 
densities found, only 417 grouped together incommensurable manifolds. The final 
line gives data for the non-arithmetic manifold with highest cusp density found. 
(The maximum possible cusp density is 0.853276..., which occurs for the figure 
eight knot complement.) 

More details of the fields occurring above are listed in the Table [3) 



degree discriminant signature minimum polynomial 



4 


1025 


0, 


2 




-x^ + 3x2 _ 


- 2x + 4 


7 


-76154488 


1, 


3 




- 2x^ + 3x5 


- 5x'^ + x^ - 8x^ - 


3 


-107 


1, 


1 


x'^ 


- x^ 3x - 


2 


6 


-1225043 


0, 


3 


x^ 


- 2x5 _ 2x^ 


+ 30x2 _ 522, ^ 29 


4 


144 


0, 


2 


x" 


-x2-Hl 






Table 


3. 


Fields 


in previous 


table. 



9. Appendix 

The indexing system used for knots and links, here and in snap, may still be 
subject to change. This is due to the difficulty of determining whether two non- 
hyperbolic links are equivalent and the consequent possibility that duplicates will 
later be discovered and removed from the tables. 
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For this reason we provide Table |4] giving the Dowker-Thistlethwaite codes of 
all the knots and links that we refer to. Since, for links, this code is not well known, 
we describe here how this works. 

Firstly there is the trivial matter of passing between the alphabetic codes used by 
snap, and their numerical forms. The Dowker code for a link with n crossings and 
k components is a permutation of the even integers 2, . . . ,2n, with possible sign 
changes, bracketed into k subsequences: e.g. (6,-8) (—10, 14, —12, —16, —2, 4). 
To express this alphabetically we encode n and k as the first two letters using a = 
1,6 = 2, etc. Then follow k letters giving the lengths of the bracketed subsequences. 
Finally there are n letters giving the sequence of even integers using a = 2, 5 = 4, 
etc. and A = —2, B = —4, etc. The alphabetic code for the above example is thus 
hbbf cDEgFHAb. 

To go from a link diagram to its Dowker code, proceed as follows. Traverse each 
component, numbering the crossings, starting with 1 on an overcrossing. When the 
first component is done, continue with consecutive numbers on the next component. 
Each crossing will receive two numbers. We can number the crossings in such a way 
that every crossing gets one even and one odd number. (This follows easily from 
the fact that we can two-colour the plane containing a link diagram.) Negate any 
even number which labels an over-crossing (for an alternating link there will not 
be any). For each odd number 1, 3, . . . , , 2?! — 1 write down the corresponding even 
number: this gives a sequence of n even numbers. A component with 2j crossings 
will have j odd numbers on it, so there will be j corresponding numbers for it in the 



code; bracket together the numbers for each component. See for example Figure 10 




Figure 10. Complete the numbering indicated on the above 
two-component link to obtain the Dowker-Thistlethwaite code 
(6,-8) (-10,14,-12,-16,-2,4). 



The reverse procedure, going from a code to a link diagram is a little bit more 
tricky, but essentially the same as the procedure for knots, described in [5]. We 
draw the first component as a knot, with extra as- yet-unconnected crossings on 
it. When adding further components we will encounter crossings with other link 
components. 
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name 


Dowker code 


name 


Dowker code 


4al 


dadbcda 


12a2937 


Icbccchjklaiefgbd 


6a5 


fbccdefacb 


12a2961 


Icbf dcf ahib ekdlgj 


6a8 


fcbbbceafbd 


12a3039 


Icbhbcfaikbj elgdh 


8a25 


hbbfcfcgahdb 


12a3169 


Iccfcdfhcj akblieg 


8a37 


hcbcccfghadeb 


12a3230 


Idbccdcfagiebkdlhj 


8a39 


hdbbbb ceagbhdf 


12a3273 


Idbddbceagbidkflhj 


8n7 


hcbcccDFgHEaB 


12a3284 


Idbfbbccaibkjldgfh 


8n8 


hcbcccdFgHEab 


12a3285 


Idbfbbcf aikbj lehdg 


8nl0 


hdbbbbcEaGBHDF 


12a3292 


Idccccdegj hkflacib 


9a73 


ibcfdfhagbice 


12a3295 


Idccccdgj akhbflice 


9a74 


ibcfdfhaibecg 


12a3296 


Idccccdgjhkbelcfia 


9a94 


icbebcdhfiabeg 


12a3300 


lebbbdbceagbidkflhj 


9n34 


icbdcccaGbHDIF 


12a3307 


Icbccbbcfaikbljdheg 


9n36 


icbebcdHfiabeG 


12a3308 


Ifbbbbbbceagbidkflhj 


10a242 


jbeefghjiaecbd 


12n706 


lalceFhGjIKlBaD 


10a277 


jcbfbceagbidjfh 


12n968 


IbbjcDFhIJLAbKEG 


10a280 


jcccddeghjibcfa 


12nll89 


IbbjceaHbKiDjgLF 


10a291 


j dbbdbceagbidj fh 


12nll90 


IbbjceaHbklDJGlf 


10a297 


jebbbbbceagbidjfh 


12nl481 


IbbjchaEGJDbkFli 


10n73 


jbbhcEFihGJAdb 


12nl848 


IbcidflaglceKBHJ 


10n74 


jbbhcfaHGIbDJE 


12n2348 


IcbbhcEaHBKiDjgLF 


10nll2 


jcbcecdFglHabJE 


12n2492 


IcbdfceaGbiDJKfLH 


10nll3 


jcbcecFaHJIBDGE 


12n2580 


IcbeecHaEGJDBkFli 


10nl30 


jcbfbceaGblDJFH 


12n2625 


IcbfdcfalJbKLEDHG 


10nl39 


jdbbcccEaHBijDgf 


12n2630 


IcbfdcfalJbKLHDEG 


10nl43 


jdbbdbcEaGBiDjfli 


12n2899 


IdbbbfcEaGBIDjkFlh 


10nl55 


jebbbbbcEaGBIDjFh 


12n2972 


IdbccdcfalJbKLEDHG 


lln318 


kbchdEfcHiJKaBG 


12n3068 


IdbddbceaGblDKFLHJ 


lln498 


kcbcfcfaHJblKEGD 


12n3078 


IdbfbbcealbKJLDGFH 


lln539 


kcbedceaHblJDGKF 


12n3088 


IdccccdeGacJBklFih 


lln545 


kcbfccdIfJabKGEH 


12n3093 


IdccccdEGHiJKlAFBc 


12al743 


Ibbjchfjialkedbg 


12n3098 


IdccccdEGhlJKLaFBC 


12a2126 


Ibcidgjclhafkbie 


12n3099 


IdccccdEGhlJKlAFBc 



Table 4. Dowkcr-Thistlethwaite codes of all knots and links men- 
tioned in this paper. 



References 

[1] C. Adams, Noncompact hyperbolic 3-orbifolds of small volume, Topology '90 (Columbus, 
OH, 1990), 1-15, Ohio State Univ. Math. Res. Inst. Publ., 1, de Gruyter, Berlin, 1992. 

[2] C. Adams, The Knot Book, W. H. Freeman and Company, New York, 1994. 

[3] I. R. Aitchison and J. H. Rubinstein, Combinatorial Cubings, Cusps and the Dodecahedral 
Knots, Topology '90 (Columbus, OH, 1990), 273-310, Ohio State Univ. Math. Res. Inst. 
Publ., 1, dc Gruyter, Berlin, 1992. 

[4] H. Akiyoshi, Finiteness of polyhedral decompositions of cusped hyperbolic manifolds ob- 
tained by the Epstein-Penner's method. Proc. Amer. Math. Soc. 129 (2001), no. 8, 2431- 
2439. 

[5] H. Akiyoshi, M. Sakuma, M. Wada and Y. Yaiuasliita, Punctured Torus Groups and 2- 
Bridge Knot Groups (I), Lecture Notes in Mathematics, Vol. 1909, Springer, 2007. 

[6] A. Borel, Commensurability classes and volumes of hyperbolic 3-manifolds. Ann. Scuola 
Norm. Sup. Pisa CI. Sci. (4) 8 (1981), no. 1, 1-33. 

[7] B. H. Bowditch, C. Maclachlan and A. W. Reid, Arithmetic hyperbolic surface bundles. 
Math. Ann. 302 (1995), no. 1, 31-60. 



32 



OLIVER GOODMAN, DAMIAN HEARD, AND CRAIG HODGSON 



[8] J. O. Button, Fibred and virtually fibred hyperbolic 3-manifolds in the censuses. Experiment. 

Math. 14 (2005), no. 2, 231-255. 
[9] P. J. Callahan, M. V. Hildebrand and J. R. Weeks, A census of cuspcd hyperbolic 3- 

manifolds. (With microfiche supplement.) Math. Comp. 68 (1999), no. 225, 321-332. 
[10] D. Coulson, O. Goodman, C. Hodgson and W. Neumann, Computing arithmetic invariants 

of 3-maiiifolds, Experimental Mathematics 9 (2000), 127-152. 
[11] D. B. A. Epstein and R. C. Penner, Euclidean decompositions of noncompact hyperbolic 

manifolds. J. Differential Geom. 27 (1988), no. 1, 67-80. 
[12] W. Floyd and A. Hatcher, Incompressible surfaces in punctured-torus bundles. Topology 

and its Applications 13 (1982), 263-282. 
[13] F. Gonzalez-Acufia and W. C. Whitten, Imbeddings of three-manifold groups. Mem. Amer. 

Math. Soc. 99 (1992), no. 474, viii-|-55 pp. 
[14] F. Gueritaud (with an appendix by D. Futer), On canonical triangulations of oncc-punctured 

torus bundles and two-bridge link complements, Geometry and Topology 10 (2006), pp. 

1239-1284. 

[15] F. Gueritaud, Geometric hyperbolique effective et triangulations ideals canoniques en di- 
mension 3, These de doctorat, Univ. Paris-sud, Orsay, 2006. 
[16] O. Goodman, Snap, the computer program, http://www.ms. unimelb.edu.au/''snap/ and 

http : // sourceforgc . net / projects / snap-pari . 
[17] O. Goodman, Tube: a computer program for studying geodesies in hyperbolie 3-manifolds, 

available from http://www.ms.unimclb.cdu.au/~snap/. 
[18] D. Heard, Orb: a computer progrurn for studying hyperbolie 3-orbifolds, available from 

http://www.ms.unimelb.edu.au/~snap/orb.html. 
[19] M. Hildebrand and J. Weeks, A computer generated census of cusped hyperbolic 3-manifolds. 

Computers and mathematics (Cambridge, MA, 1989), 53-59, Springer, New York, 1989. 
[20] J. Hoste and M. Thistlethwaite, Knotseape, http://www.math.utk.edu/~morwen/knotscape.html. 
[21] M. Lackcnby, The canonical decomposition of once-punctured torus bundles, Comment. 

Math. Helv. 78 (2003) 363-384. 
[22] A. M. Macbcath, Commcnsurability of co-compact three-dimensional hyperbolic groups. 

Duke Math. J. 50 (1983), no. 4, 1245-1253. Erratum: Duke Math. J. 56 (1988), no. 1, 219. 
[23] C. Maclachlan and A. Reid, The arithmetic of hyperbolic 3-manifolds, Springer- Verlag, New 

York, 2003. 

[24] G. A. Margulis, Discrete subgroups of semisimple Lie groups, Ergebnisse der Mathematik 

und ihrer Grenzgebiete 17, Springer- Verlag, Berlin, 1991. 
[25] T. H. Marshall and G. J. Martin, Minimal Co- Volume Hyperbolic Lattices H, in preparation. 
[26] R. Meyerhoff, The cusped hyperbolic 3-orbifold of minimal volume, Bull. Amer. Math. Soc. 

13 (1985), 154-156. 

[27] W. D. Neumann and A. W. I^cid. Arithmetic of Hyperbolic Manifolds, Topology '90 (Colum- 
bus, OH, 1990), 273-310, Ohio State Univ. Math. Res. Inst. Publ., 1, de Gruytcr, Berlin, 
1992. 

[28] W. D. Neumann and A. W. Reid, Notes on Adams' small volume orbifolds. Topology '90 
(Columbus, OH, 1990), 311-314, Ohio State Univ. Math. Res. Inst. Publ., 1, de Gruyter, 
BerUn, 1992. 

[29] A. W. Reid, A note on trace-fields of Kleinian groups. Bull. London Math. Soc. 22 (1990), 
no. 4, 349-352. 

[30] M. Sakuma and J. Weeks, Examples of canonical decompositions of hyperbolic link comple- 
ments. Japan. J. Math. (N.S.) 21 (1995), no. 2, 393-439. 

[31] M. Sakuma and J. Weeks, The generalized tilt formula. Geom. Dcdicata 55 (1995), no. 2, 
115-123. 

[32] W. P. Thurston, The geometry and topology of three-manifolds, Princeton University Math. 

Dept. (1978). Also available at http://msri.org/publications/books/gt3m. 
[33] J. R. Weeks, Convex hulls and isometries of cusped hyperbolic 3-manifolds. Topology Appl. 

52 (1993), no. 2, 127-149. 
[34] J. Weeks, SnapPea, the computer program, available from 

http: / /geometrygames.org/SnapPea/index.html. 
[35] R. Zimmer, Ergodic theory and semi-simple Lie groups, Birkhauser, Boston, 1984. 

Department of Mathematics and Statistics, University of Melbourne, Parkville, Vic- 
toria 3010, Australia 

E-mail address: oagSoptusnet.com.au, damian.heardagmail . com, cdhOms .unimelb . edu. au 



